Recurrent SARS-CoV-2 mutations at Spike D796 evade antibodies from pre-Omicron convalescent and vaccinated subjects.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lineages of the Omicron variant rapidly became dominant in early 2022 and frequently cause human infections despite vaccination or prior infection with other variants. In addition to antibody-evading mutations in the receptor-binding domain, Omicron features amino acid mutations elsewhere in the Spike protein; however, their effects generally remain ill defined. The Spike D796Y substitution is present in all Omicron sub-variants and occurs at the same site as a mutation (D796H) selected during viral evolution in a chronically infected patient. Here, we map antibody reactivity to a linear epitope in the Spike protein overlapping position 796. We show that antibodies binding this region arise in pre-Omicron SARS-CoV-2 convalescent and vaccinated subjects but that both D796Y and D796H abrogate their binding. These results suggest that D796Y contributes to the fitness of Omicron in hosts with pre-existing immunity to other variants of SARS-CoV-2 by evading antibodies targeting this site.IMPORTANCESevere acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has evolved substantially through the coronavirus disease 2019 (COVID-19) pandemic: understanding the drivers and consequences of this evolution is essential for projecting the course of the pandemic and developing new countermeasures. Here, we study the immunological effects of a particular mutation present in the Spike protein of all Omicron strains and find that it prevents the efficient binding of a class of antibodies raised by pre-Omicron vaccination and infection. These findings reveal a novel consequence of a poorly understood Omicron mutation and shed light on the drivers and effects of SARS-CoV-2 evolution.

Global phylogenomic diversity of Brucella abortus: spread of a dominant lineage.

Brucella abortus is a globally important zoonotic pathogen largely found in cattle hosts and is typically transmitted to humans through contaminated dairy products or contact with diseased animals. Despite the long, shared history of cattle and humans, little is known about how trade in cattle has spread this pathogen throughout the world. Whole genome sequencing provides unparalleled resolution to investigate the global evolutionary history of a bacterium such as B. abortus by providing phylogenetic resolution that has been unobtainable using other methods. We report on large-scale genome sequencing and analysis of B. abortus collected globally from cattle and 16 other hosts from 52 countries. We used single nucleotide polymorphisms (SNPs) to identify genetic variation in 1,074 B. abortus genomes and using maximum parsimony generated a phylogeny that identified four major clades. Two of these clades, clade A (median date 972 CE; 95% HPD, 781-1142 CE) and clade B (median date 150 BCE; 95% HPD, 515 BCE-164 CE), were exceptionally diverse for this species and are exclusively of African origin where provenance is known. The third clade, clade C (median date 949 CE; 95% HPD, 766-1102 CE), had most isolates coming from a broad swath of the Middle East, Europe, and Asia, also had relatively high diversity. Finally, the fourth major clade, clade D (median date 1467 CE; 95% HPD, 1367-1553 CE) comprises the large majority of genomes in a dominant but relatively monomorphic group that predominantly infects cattle in Europe and the Americas. These data are consistent with an African origin for B. abortus and a subsequent spread to the Middle East, Europe, and Asia, probably through the movement of infected cattle. We hypothesize that European arrival to the Americas starting in the 15th century introduced B. abortus from Western Europe through the introduction of a few common cattle breeds infected with strains from clade D. These data provide the foundation of a comprehensive global phylogeny of this important zoonotic pathogen that should be an important resource in human and veterinary epidemiology.

Towards a post-pandemic future for global pathogen genome sequencing.

Pathogen genome sequencing has become a routine part of our response to active outbreaks of infectious disease and should be an important part of our preparations for future epidemics. In this Essay, we discuss the innovations that have enabled routine pathogen genome sequencing, as well as how genome sequences can be used to understand and control the spread of infectious disease. We also explore the impact of the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) pandemic on the field of pathogen genomics and outline the challenges we must address to further improve the utility of pathogen genome sequencing in the future.

Virome-wide detection of natural infection events and the associated antibody dynamics using longitudinal highly-multiplexed serology.

Current methods for detecting infections either require a sample collected from an actively infected site, are limited in the number of agents they can query, and/or yield no information on the immune response. Here we present an approach that uses temporally coordinated changes in highly-multiplexed antibody measurements from longitudinal blood samples to monitor infection events at sub-species resolution across the human virome. In a longitudinally-sampled cohort of South African adolescents representing >100 person-years, we identify >650 events across 48 virus species and observe strong epidemic effects, including high-incidence waves of Aichivirus A and the D68 subtype of Enterovirus D earlier than their widespread circulation was appreciated. In separate cohorts of adults who were sampled at higher frequency using self-collected dried blood spots, we show that such events temporally correlate with symptoms and transient inflammatory biomarker elevations, and observe the responding antibodies to persist for periods ranging from ≤1 week to >5 years. Our approach generates a rich view of viral/host dynamics, supporting novel studies in immunology and epidemiology.

PepSeq: a fully in vitro platform for highly multiplexed serology using customizable DNA-barcoded peptide libraries.

PepSeq is an in vitro platform for building and conducting highly multiplexed proteomic assays against customizable targets by using DNA-barcoded peptides. Starting with a pool of DNA oligonucleotides encoding peptides of interest, this protocol outlines a fully in vitro and massively parallel procedure for synthesizing the encoded peptides and covalently linking each to a corresponding cDNA tag. The resulting libraries of peptide/DNA conjugates can be used for highly multiplexed assays that leverage high-throughput sequencing to profile the binding or enzymatic specificities of proteins of interest. Here, we describe the implementation of PepSeq for fast and cost-effective epitope-level analysis of antibody reactivity across hundreds of thousands of peptides from <1 µl of serum or plasma input. This protocol includes the design of the DNA oligonucleotide library, synthesis of DNA-barcoded peptide constructs, binding of constructs to sample, preparation for sequencing and data analysis. Implemented in this way, PepSeq can be used for a number of applications, including fine-scale mapping of antibody epitopes and determining a subject's pathogen exposure history. The protocol is divided into two main sections: (i) design and synthesis of DNA-barcoded peptide libraries and (ii) use of libraries for highly multiplexed serology. Once oligonucleotide templates are in hand, library synthesis takes 1-2 weeks and can provide enough material for hundreds to thousands of assays. Serological assays can be conducted in 96-well plates and generate sequencing data within a further ~4 d. A suite of software tools, including the PepSIRF package, are made available to facilitate the design of PepSeq libraries and analysis of assay data.

Highly Multiplexed Serology for Nonhuman Mammals.

Emerging infectious diseases represent a serious and ongoing threat to humans. Most emerging viruses are maintained in stable relationships with other species of animals, and their emergence within the human population results from cross-species transmission. Therefore, if we want to be prepared for the next emerging virus, we need to broadly characterize the diversity and ecology of viruses currently infecting other animals (i.e., the animal virosphere). High-throughput metagenomic sequencing has accelerated the pace of virus discovery. However, molecular assays can detect only active infections and only if virus is present within the sampled fluid or tissue at the time of collection. In contrast, serological assays measure long-lived antibody responses to infections, which can be detected within the blood, regardless of the infected tissues. Therefore, serological assays can provide a complementary approach for understanding the circulation of viruses, and while serological assays have historically been limited in scope, recent advancements allow thousands to hundreds of thousands of antigens to be assessed simultaneously using <1 µL of blood (i.e., highly multiplexed serology). The application of highly multiplexed serology for the characterization of the animal virosphere is dependent on the availability of reagents that can be used to capture or label antibodies of interest. Here, we evaluate the utility of commercial immunoglobulin-binding proteins (protein A and protein G) to enable highly multiplexed serology in 25 species of nonhuman mammals, and we describe a competitive fluorescence-linked immunosorbent assay (FLISA) that can be used as an initial screen for choosing the most appropriate capture protein for a given host species. IMPORTANCE Antibodies are generated in response to infections with viruses and other pathogens, and they help protect against future exposures. Mature antibodies are long lived, are highly specific, and can bind to their protein targets with high affinity. Thus, antibodies can also provide information about an individual's history of viral exposures, which has important applications for understanding the epidemiology and etiology of disease. In recent years, there have been large advances in the available methods for broadly characterizing antibody-binding profiles, but thus far, these have been utilized primarily with human samples only. Here, we demonstrate that commercial antibody-binding reagents can facilitate modern antibody assays for a wide variety of mammalian species, and we describe an inexpensive and fast approach for choosing the best reagent for each animal species. By studying antibody-binding profiles in captive and wild animals, we can better understand the distribution and prevalence of viruses that could spill over into humans.

COVID-19 vaccination elicits an evolving, cross-reactive antibody response to epitopes conserved with endemic coronavirus spike proteins.

The COVID-19 pandemic has triggered the first widespread vaccination campaign against a coronavirus. Many vaccinated subjects are previously naive to SARS-CoV-2; however, almost all have previously encountered other coronaviruses (CoVs), and the role of this immunity in shaping the vaccine response remains uncharacterized. Here, we use longitudinal samples and highly multiplexed serology to identify mRNA-1273 vaccine-induced antibody responses against a range of CoV Spike epitopes, in both phylogenetically conserved and non-conserved regions. Whereas reactivity to SARS-CoV-2 epitopes shows a delayed but progressive increase following vaccination, we observe distinct kinetics for the endemic CoV homologs at conserved sites in Spike S2: these become detectable sooner and decay at later time points. Using homolog-specific antibody depletion and alanine-substitution experiments, we show that these distinct trajectories reflect an evolving cross-reactive response that can distinguish rare, polymorphic residues within these epitopes. Our results reveal mechanisms for the formation of antibodies with broad reactivity against CoVs.

The population genetics of the causative agent of snake fungal disease indicate recent introductions to the USA.

Snake fungal disease (SFD; ophidiomycosis), caused by the pathogen Ophidiomyces ophiodiicola (Oo), has been documented in wild snakes in North America and Eurasia, and is considered an emerging disease in the eastern United States of America. However, a lack of historical disease data has made it challenging to determine whether Oo is a recent arrival to the USA or whether SFD emergence is due to other factors. Here, we examined the genomes of 82 Oo strains to determine the pathogen's history in the eastern USA. Oo strains from the USA formed a clade (Clade II) distinct from European strains (Clade I), and molecular dating indicated that these clades diverged too recently (approximately 2,000 years ago) for transcontinental dispersal of Oo to have occurred via natural snake movements across Beringia. A lack of nonrecombinant intermediates between clonal lineages in Clade II indicates that Oo has actually been introduced multiple times to North America from an unsampled source population, and molecular dating indicates that several of these introductions occurred within the last few hundred years. Molecular dating also indicated that the most common Clade II clonal lineages have expanded recently in the USA, with time of most recent common ancestor mean estimates ranging from 1985 to 2007 CE. The presence of Clade II in captive snakes worldwide demonstrates a potential mechanism of introduction and highlights that additional incursions are likely unless action is taken to reduce the risk of pathogen translocation and spillover into wild snake populations.

COVID-19 vaccination recruits and matures cross-reactive antibodies to conserved epitopes in endemic coronavirus Spike proteins.

The COVID-19 pandemic has triggered the first widespread vaccination campaign against a coronavirus. Most vaccinated subjects are naïve to SARS-CoV-2, however almost all have previously encountered other coronaviruses (CoVs) and the role of this immunity in shaping the vaccine response remains uncharacterized. Here we use longitudinal samples and highly-multiplexed serology to identify mRNA-1273 vaccine-induced antibody responses against a range of CoV Spike epitopes and in both phylogenetically conserved and non-conserved regions. Whereas reactivity to SARS-CoV-2 epitopes showed a delayed but progressive increase following vaccination, we observed distinct kinetics for the endemic CoV homologs at two conserved sites in Spike S2: these became detectable sooner, and decayed at later timepoints. Using homolog-specific depletion and alanine-substitution experiments, we show that these distinctly-evolving specificities result from cross-reactive antibodies as they mature against rare, polymorphic residues within these epitopes. Our results reveal mechanisms for the formation of antibodies with broad reactivity against CoVs.

Comparative Genomics Analyses Support the Reclassification of Bisgaard Taxon 40 as Mergibacter gen. nov., With Mergibacter septicus sp. nov. as Type Species: Novel Insights Into the Phylogeny and Virulence Factors of a Pasteurellaceae Family Member Associated With Mortality Events in Seabirds.

The Pasteurellaceae family has been associated with fatal diseases in numerous avian species. Several new taxa within this family, including Bisgaard taxon 40, have been recently described in wild birds, but their genomic characteristics and pathogenicity are not well understood. We isolated Bisgaard taxon 40 from four species of seabirds, including one sampled during a mass, multi-species mortality event in Florida, United States. Here, we present a comprehensive phenotypic and genetic characterization of Bisgaard taxon 40 and comparative genomic analysis with reference strains from the Pasteurellaceae family, aiming at determining its phylogenetic position, antimicrobial susceptibility profile, and identifying putative virulence factors. In silico multilocus sequence-based and whole-genome-based phylogenetic analysis clustered all Bisgaard taxon 40 strains together on a distinct branch separated from the other members of the Pasteurellaceae family, indicating that Bisgaard taxon 40 could represent a new genus. These findings were further supported by protein similarity analyses using the concatenation of 31 conserved proteins and other taxonomic approaches such as the percentage of conserved protein test. Additionally, several putative virulence factors were identified, including those associated with adhesion (capsule, ompA, ompH) and colonization (exbD, fur, galU, galE, lpxA, lpxC, and kdsA) of the host and a cytolethal distending toxin (cdt), which may have played a role in disease development leading to the mortality event. Considerably low minimum inhibitory concentrations (MICs) were found for all the drugs tested, in concordance with the absence of antimicrobial resistance genes in these genomes. The novel findings of this study highlight genomic and phenotypic characteristics of this bacterium, providing insights into genome evolution and pathogenicity. We propose a reclassification of these organisms within the Pasteurellaceae family, designated as Mergibacter gen. nov., with Mergibacter septicus sp. nov. as the type species. The type strain is Mergibacter septicus A25201T (=DSM 112696).

Genomic signatures for predicting the zoonotic potential of novel viruses.

Powered by metagenomics, viral discovery is outpacing our capacity for the downstream characterization needed to fully assess zoonotic potential. A study published in PLOS Biology uses machine learning to prioritize novel viruses based only on genomic signatures.

Multiplexed Simian Immunodeficiency Virus-Specific Paired RNA-Guided Cas9 Nickases Inactivate Proviral DNA.

Human and simian immunodeficiency virus (HIV and SIV) infections establish lifelong reservoirs of cells harboring an integrated proviral genome. Genome editing CRISPR-associated Cas9 nucleases, combined with SIV-specific guiding RNA (gRNA) molecules, inactivate integrated provirus DNA in vitro and in animal models. We generated RNA-guided Cas9 nucleases (RGNu) and nickases (RGNi) targeting conserved SIV regions with no homology in the human or rhesus macaque genome. Assays in cells cotransfected with SIV provirus and plasmids coding for RGNus identified SIV long terminal repeat (LTR), trans-activation response (TAR) element, and ribosome slip site (RSS) regions as the most effective at virus suppression; RGNi targeting these regions inhibited virus production significantly. Multiplex plasmids that coexpressed these three RGNu (Nu3), or six (three pairs) RGNi (Ni6), were more efficient at virus suppression than any combination of individual RGNu and RGNi plasmids. Both Nu3 and Ni6 plasmids were tested in lymphoid cells chronically infected with SIVmac239, and whole-genome sequencing was used to determine on- and off-target mutations. Treatment with these all-in-one plasmids resulted in similar levels of mutations of viral sequences from the cellular genome; Nu3 induced indels at the 3 SIV-specific sites, whereas for Ni6 indels were present at the LTR and TAR sites. Levels of off-target effects detected by two different algorithms were indistinguishable from background mutations. In summary, we demonstrate that Cas9 nickase in association with gRNA pairs can specifically eliminate parts of the integrated provirus DNA; also, we show that careful design of an all-in-one plasmid coding for 3 gRNAs and Cas9 nuclease inhibits SIV production with undetectable off-target mutations, making these tools a desirable prospect for moving into animal studies. IMPORTANCE Our approach to HIV cure, utilizing the translatable SIV/rhesus macaque model system, aims at provirus inactivation and its removal with the least possible off-target side effects. We developed single molecules that delivered either three truncated SIV-specific gRNAs along with Cas9 nuclease or three pairs of SIV-specific gRNAs (six individual gRNAs) along with Cas9 nickase to enhance efficacy of on-target mutagenesis. Whole-genome sequencing demonstrated effective SIV sequence mutation and inactivation and the absence of demonstrable off-target mutations. These results open the possibility to employ Cas9 variants that introduce single-strand DNA breaks to eliminate integrated proviral DNA.

Inference of Nipah virus evolution, 1999-2015.

Despite near-annual human outbreaks of Nipah virus (NiV) disease in Bangladesh, typically due to individual spillover events from the local bat population, only twenty whole-genome NiV sequences exist from humans and ten from bats. NiV whole-genome sequences from annual outbreaks have been challenging to generate, primarily due to the low viral load in human throat swab and serum specimens. Here, we used targeted enrichment with custom NiV-specific probes and generated thirty-five additional unique full-length genomic sequences directly from human specimens and viral isolates. We inferred the temporal and geographic evolutionary history of NiV in Bangladesh and expanded a tool to visualize NiV spatio-temporal spread from a Bayesian continuous diffusion analysis. We observed that strains from Bangladesh segregated into two distinct clades that have intermingled geographically in Bangladesh over time and space. As these clades expanded geographically and temporally, we did not observe evidence for significant branch and site-specific selection, except for a single site in the Henipavirus L polymerase. However, the Bangladesh 1 and 2 clades are differentiated by mutations initially occurring in the polymerase, with additional mutations accumulating in the N, G, F, P, and L genes on external branches. Modeling the historic geographical and temporal spread demonstrates that while widespread, NiV does not exhibit significant genetic variation in Bangladesh. Thus, future public health measures should address whether NiV within in the bat population also exhibits comparable genetic variation, if zoonotic transmission results in a genetic bottleneck and if surveillance techniques are detecting only a subset of NiV.

Epitope-resolved profiling of the SARS-CoV-2 antibody response identifies cross-reactivity with endemic human coronaviruses.

The SARS-CoV-2 proteome shares regions of conservation with endemic human coronaviruses (CoVs), but it remains unknown to what extent these may be cross-recognized by the antibody response. Here, we study cross-reactivity using a highly multiplexed peptide assay (PepSeq) to generate an epitope-resolved view of IgG reactivity across all human CoVs in both COVID-19 convalescent and negative donors. PepSeq resolves epitopes across the SARS-CoV-2 Spike and Nucleocapsid proteins that are commonly targeted in convalescent donors, including several sites also recognized in some uninfected controls. By comparing patterns of homologous reactivity between CoVs and using targeted antibody-depletion experiments, we demonstrate that SARS-CoV-2 elicits antibodies that cross-recognize pandemic and endemic CoV antigens at two Spike S2 subunit epitopes. We further show that these cross-reactive antibodies preferentially bind endemic homologs. Our findings highlight sites at which the SARS-CoV-2 response appears to be shaped by previous CoV exposures and which have the potential to raise broadly neutralizing responses.

Evolution of Antibiotic Resistance in Surrogates of Francisella tularensis (LVS and Francisella novicida): Effects on Biofilm Formation and Fitness.

Francisella tularensis, the causative agent of tularemia, is capable of causing disease in a multitude of mammals and remains a formidable human pathogen due to a high morbidity, low infectious dose, lack of a FDA approved vaccine, and ease of aerosolization. For these reasons, there is concern over the use of F. tularensis as a biological weapon, and, therefore, it has been classified as a Tier 1 select agent. Fluoroquinolones and aminoglycosides often serve as the first line of defense for treatment of tularemia. However, high levels of resistance to these antibiotics has been observed in gram-negative bacteria in recent years, and naturally derived resistant Francisella strains have been described in the literature. The acquisition of antibiotic resistance, either natural or engineered, presents a challenge for the development of medical countermeasures. In this study, we generated a surrogate panel of antibiotic resistant F. novicida and Live Vaccine Strain (LVS) by selection in the presence of antibiotics and characterized their growth, biofilm capacity, and fitness. These experiments were carried out in an effort to (1) assess the fitness of resistant strains; and (2) identify new targets to investigate for the development of vaccines or therapeutics. All strains exhibited a high level of resistance to either ciprofloxacin or streptomycin, a fluoroquinolone and aminoglycoside, respectively. Whole genome sequencing of this panel revealed both on-pathway and off-pathway mutations, with more mutations arising in LVS. For F. novicida, we observed decreased biofilm formation for all ciprofloxacin resistant strains compared to wild-type, while streptomycin resistant isolates were unaffected in biofilm capacity. The fitness of representative antibiotic resistant strains was assessed in vitro in murine macrophage-like cell lines, and also in vivo in a murine model of pneumonic infection. These experiments revealed that mutations obtained by these methods led to nearly all ciprofloxacin resistant Francisella strains tested being completely attenuated while mild attenuation was observed in streptomycin resistant strains. This study is one of the few to examine the link between acquired antibiotic resistance and fitness in Francisella spp., as well as enable the discovery of new targets for medical countermeasure development.

An Early Pandemic Analysis of SARS-CoV-2 Population Structure and Dynamics in Arizona.

In December of 2019, a novel coronavirus, SARS-CoV-2, emerged in the city of Wuhan, China, causing severe morbidity and mortality. Since then, the virus has swept across the globe, causing millions of confirmed infections and hundreds of thousands of deaths. To better understand the nature of the pandemic and the introduction and spread of the virus in Arizona, we sequenced viral genomes from clinical samples tested at the TGen North Clinical Laboratory, the Arizona Department of Health Services, and those collected as part of community surveillance projects at Arizona State University and the University of Arizona. Phylogenetic analysis of 84 genomes from across Arizona revealed a minimum of 11 distinct introductions inferred to have occurred during February and March. We show that >80% of our sequences descend from strains that were initially circulating widely in Europe but have since dominated the outbreak in the United States. In addition, we show that the first reported case of community transmission in Arizona descended from the Washington state outbreak that was discovered in late February. Notably, none of the observed transmission clusters are epidemiologically linked to the original travel-related case in the state, suggesting successful early isolation and quarantine. Finally, we use molecular clock analyses to demonstrate a lack of identifiable, widespread cryptic transmission in Arizona prior to the middle of February 2020.IMPORTANCE As the COVID-19 pandemic swept across the United States, there was great differential impact on local and regional communities. One of the earliest and hardest hit regions was in New York, while at the same time Arizona (for example) had low incidence. That situation has changed dramatically, with Arizona now having the highest rate of disease increase in the country. Understanding the roots of the pandemic during the initial months is essential as the pandemic continues and reaches new heights. Genomic analysis and phylogenetic modeling of SARS-COV-2 in Arizona can help to reconstruct population composition and predict the earliest undetected introductions. This foundational work represents the basis for future analysis and understanding as the pandemic continues.

Epitope-resolved profiling of the SARS-CoV-2 antibody response identifies cross-reactivity with an endemic human CoV.

A high-resolution understanding of the antibody response to SARS-CoV-2 is important for the design of effective diagnostics, vaccines and therapeutics. However, SARS-CoV-2 antibody epitopes remain largely uncharacterized, and it is unknown whether and how the response may cross-react with related viruses. Here, we use a multiplexed peptide assay ('PepSeq') to generate an epitope-resolved view of reactivity across all human coronaviruses. PepSeq accurately detects SARS-CoV-2 exposure and resolves epitopes across the Spike and Nucleocapsid proteins. Two of these represent recurrent reactivities to conserved, functionally-important sites in the Spike S2 subunit, regions that we show are also targeted for the endemic coronaviruses in pre-pandemic controls. At one of these sites, we demonstrate that the SARS-CoV-2 response strongly and recurrently cross-reacts with the endemic virus hCoV-OC43. Our analyses reveal new diagnostic and therapeutic targets, including a site at which SARS-CoV-2 may recruit common pre-existing antibodies and with the potential for broadly-neutralizing responses.

Persistence of Brucella abortus lineages revealed by genomic characterization and phylodynamic analysis.

Brucellosis, caused by Brucella abortus, is a major disease of cattle and humans worldwide distributed. Eradication and control of the disease has been difficult in Central and South America, Central Asia, the Mediterranean and the Middle East. Epidemiological strategies combined with phylogenetic methods provide the high-resolution power needed to study relationships between surveillance data and pathogen population dynamics, using genetic diversity and spatiotemporal distributions. This information is crucial for prevention and control of disease spreading at a local and worldwide level. In Costa Rica (CR), the disease was first reported at the beginning of the 20th century and has not been controlled despite many efforts. We characterized 188 B. abortus isolates from CR recovered from cattle, humans and water buffalo, from 2003 to 2018, and whole genome sequencing (WGS) was performed in 95 of them. They were also assessed based on geographic origin, date of introduction, and phylogenetic associations in a worldwide and national context. Our results show circulation of five B. abortus lineages (I to V) in CR, phylogenetically related to isolates from the United States, United Kingdom, and South America. Lineage I was dominant and probably introduced at the end of the 19th century. Lineage II, represented by a single isolate from a water buffalo, clustered with a Colombian sample, and was likely introduced after 1845. Lineages III and IV were likely introduced during the early 2000s. Fourteen isolates from humans were found within the same lineage (lineage I) regardless of their geographic origin within the country. The main CR lineages, introduced more than 100 years ago, are widely spread throughout the country, in contrast to new introductions that seemed to be more geographically restricted. Following the brucellosis prevalence and the farming practices of several middle- and low-income countries, similar scenarios could be found in other regions worldwide.

Reproducibly sampling SARS-CoV-2 genomes across time, geography, and viral diversity.

The COVID-19 pandemic has led to a rapid accumulation of SARS-CoV-2 genomes, enabling genomic epidemiology on local and global scales. Collections of genomes from resources such as GISAID must be subsampled to enable computationally feasible phylogenetic and other analyses. We present genome-sampler, a software package that supports sampling collections of viral genomes across multiple axes including time of genome isolation, location of genome isolation, and viral diversity. The software is modular in design so that these or future sampling approaches can be applied independently and combined (or replaced with a random sampling approach) to facilitate custom workflows and benchmarking. genome-sampler is written as a QIIME 2 plugin, ensuring that its application is fully reproducible through QIIME 2's unique retrospective data provenance tracking system. genome-sampler can be installed in a conda environment on macOS or Linux systems. A complete default pipeline is available through a Snakemake workflow, so subsampling can be achieved using a single command. genome-sampler is open source, free for all to use, and available at https://caporasolab.us/genome-sampler. We hope that this will facilitate SARS-CoV-2 research and support evaluation of viral genome sampling approaches for genomic epidemiology.

Lassa virus circulating in Liberia: a retrospective genomic characterisation.

BACKGROUND: An alarming rise in reported Lassa fever cases continues in west Africa. Liberia has the largest reported per capita incidence of Lassa fever cases in the region, but genomic information on the circulating strains is scarce. The aim of this study was to substantially increase the available pool of data to help foster the generation of targeted diagnostics and therapeutics. METHODS: Clinical serum samples collected from 17 positive Lassa fever cases originating from Liberia (16 cases) and Guinea (one case) within the past decade were processed at the Liberian Institute for Biomedical Research using a targeted-enrichment sequencing approach, producing 17 near-complete genomes. An additional 17 Lassa virus sequences (two from Guinea, seven from Liberia, four from Nigeria, and four from Sierra Leone) were generated from viral stocks at the US Centers for Disease Control and Prevention (Atlanta, GA) from samples originating from the Mano River Union (Guinea, Liberia, and Sierra Leone) region and Nigeria. Sequences were compared with existing Lassa virus genomes and published Lassa virus assays. FINDINGS: The 23 new Liberian Lassa virus genomes grouped within two clades (IV.A and IV.B) and were genetically divergent from those circulating elsewhere in west Africa. A time-calibrated phylogeographic analysis incorporating the new genomes suggests Liberia was the entry point of Lassa virus into the Mano River Union region and estimates the introduction to have occurred between 300-350 years ago. A high level of diversity exists between the Liberian Lassa virus genomes. Nucleotide percent difference between Liberian Lassa virus genomes ranged up to 27% in the L segment and 18% in the S segment. The commonly used Lassa Josiah-MGB assay was up to 25% divergent across the target sites when aligned to the Liberian Lassa virus genomes. INTERPRETATION: The large amount of novel genomic diversity of Lassa virus observed in the Liberian cases emphasises the need to match deployed diagnostic capabilities with locally circulating strains and underscores the importance of evaluating cross-lineage protection in the development of vaccines and therapeutics. FUNDING: Defense Biological Product Assurance Office of the US Department of Defense and the Armed Forces Health Surveillance Branch and its Global Emerging Infections Surveillance and Response Section.